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Abstract 

The 0{n) Gross-Neveu model for n < 2 presents a massless phase that can be 
characterized by right-left mover scattering processes. The limit n — > describes 
the on-shell properties of the random bond Ising model. 
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1. Aim of this paper is to discuss the S'-matrix formulation of the 0(n) Gross- Neveu (GN) 
model for n < 2, and in particular the scattering theory associated to the limit n — > 0. 
There is a well-defined physical problem related to this limit, which is the analysis of the 
influence of random impurities, induced by thermal fluctuations, on the critical behaviour 
of pure homogeneous systems. It is therefore useful to briefly remind the relationship 
between the GN model and the random systems. 

For "random temperature" kinds of impurities, i.e. defects or dislocations in the 
material induced by thermal fluctuations, there is a simple criterion Jj]] for extimating the 
relevance of weak disorder on a critical system. According to this criterion, the effect of 
the disorder depends on the sign of the specific heat critical index a of the pure material. 
For a > the impurities are expected to completely suppress the long range fluctuations 
of the pure system, canceling all singularities in the thermodynamical quantities. On the 
contrary, for a < 0, the impurities may produce a shift of the critical temperature but 
they do not affect the critical behaviour, i.e. the critical exponents are the same as in 
the pure system. The marginal case a = is special and must be separately analysed. 
This situation occurs in the two-dimensional Ising model with random bond distribution, 
and it has been initially considered by Dotsenko and Dotsenko ||. In particular, these 
authors have shown that near the critical temperature, the class of universality of the 
random bond Ising model is described by the 0(n) Gross- Neveu model in the limit n — ► 
(see also || f|, [5|). The mapping is realized as follows. The two-dimensional homogeneous 
Ising model can be described in the continuum limit by a massive Majorana fermion 
where the mass m is a linear measurament of the deviation of the temperature T from 
the critical value T c . The partition function for the pure system is given by 



Suppose now that the temperature, instead of being uniform on the whole lattice, is 
allowed to vary from point to point but staying close, in average, to the critical value. 
The mass parameter becomes then a random variable m(x), with a probability distribution 
P[m(x)} which is supposed to be gaussian, P[m(x)} ~ exp [— m 2 (x)/2g 2 ]. Assuming that 
the time scale of the spin-flips is much faster than the typical time scale in which the local 
temperatures are updated, one is led to consider the quenched average on the disorder 
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variables, i.e. the average of the free energy rather than the average of the partition 
function 

hZrei = /lnZ[m(x)]P[m(x)]^) . (2) 

A standard method to get around the difficulty presented by this computation is to use 
the replica trick ||, which consists in replacing 

In Z — > lim 

in the integral of the right hand side of (0). In this way, taking the quenched average 
for the random bond Ising model is equivalent to solve the dynamics of the Gross-Neveu 
model in the limit n — > 

Ih~Z~ J[d&i]exp[-S\ , (3) 



where 



S = d 2 x 



i = l,...,n. (4) 



2. It is well known that the Gross-Neveu model |7] presents different properties for n > 2 
and for n < 2 ( for n = 2 it coincides with the massless Thirring model, which is a 
Conformal Field Theory). This can be seen, for instance, by considering the j3 function 
of the theory, that at the first order reads || 

m Jm 9{m) = mM)) = ~ {n ~ 2) C + ' ' ' (5) 

(M is a mass scale). For n > 2 the theory is asymptotically free in the ultraviolet limit and 
there is a dynamical mass generation. The chirality symmetry of the original lagrangian 
(fD is spontaneously broken and the fermionic fields are actually massive. The resulting 
theory is integrable and the exact elastic S-matrix has been computed in [[], In 



addition to the elementary fermions appearing in the lagrangian (|J), the final spectrum 
of the theory comes out to have quite a rich structure of bound states. 

The situation is rather different when n < 2. Since the (3 function is now positive 
around the origin, the GN model becomes an infrared asymptotically free theory and 
therefore a dynamical mass generation for the physical particles cannot occur . In this 
case, we can extract the spectrum by directly looking at the infrared limit of the lagrangian 



(|]): it simply consists in n copies of massless Majorana fermions and no additional low- 
energy excitations are expected. Since the /3-function is not zero, the model is massless but 
not conformally invariant and therefore the mass scale M entering the Renormalization 
Group equation (|5|) characterises the cross-over of the theory going from the infrared to the 
ultraviolet region. This massless phase of the GN model consists then of n right-moving 
particles R a (p) with the energy spectrum E = p (p > 0) and n left- moving particles L a (p) 
(p < 0) with E = —p. These dispersion relations can be conveniently parameterized in 
terms of the mass scale M and the rapidity variable 9 as follows: E = p = ^e e for the 
right-movers and E = —p = ^e~ e for the left-movers. 

3. Assuming that the integrability of the GN model also holds for n < 2, we can compute 
the exact elastic S-matrix associated to this massless phase of the model (for the definition 
of massless scattering and its properties, we refer the reader to the original references 



n, P2"|). In the scattering processes of integrable massless theories, the numbers of left 
and right particles are separately conserved and there is a factorization of the amplitudes. 
Hence we can restrict our attention only to the two-body processes. For the infrared 
asymptotic freedom of the model, we expect no scattering of the left-left or right-right 
particles. Therefore, in these two sectors the S'-matrix is simply —1. On the contrary, in 
the right-left sector, the S'-matrix is defined by the commutation relation 

R a (9i)L b (6 2 ) = S c a t(9i - 6 2 )L d {9 2 )Rc(9i) • (6) 

The S- matrix has to respect the 0(n) invariance of the Lagrangian (|J) and can conse- 
quently be decomposed as 

S c a f(0) = Sa b 5 cd cn(Q) + S c a Sf a 2 (6) + 5 d a 5 c b a 3 (6) . (7) 
Equivalently, we may consider the S-matrix for the channels of definite isospin, 

^isoscalar 

iVcxi + a 2 + cr 3 ; 

Cantisym = <?2 — °3 i (8) 
O'sym = CT2 + 0"3 • 
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The above amplitudes have to satisfy the unitarity and crossing symmetry equations, 
which for the massless scattering are given respectively by 

a 2 (0)a*(0) + a 3 (0)a*(0) = l ; 
n <7 X (0K(0) + <ri(0K(0) + a x (0)a*(0) + a 2 {9)a{{9) + a z {9)a{{9) = ; (9) 
a 2 (9)a* 3 (9) + a 3 (9)a* 2 (9) = 0, 

and 

a* x {9) =a,(m + 9) ; 
^1(0) =(J 2 (i7r + 0) . 

The associativity condition for the algebra ([]) gives rise to the Yang-Baxter equations, 
whose solution is given by 

with the parameter A fixed by the crossing symmetry equations (0) to be 

A = — • (12) 
n - 2 v 7 

Inserting (|TTD into eq. (|9]) and using (0), we arrive to the combined unitarity-crossing 
equation satisfied by <r 2 , 

a 2 (9)a 2 (i7r + 9) = . (13) 

The minimal solution is given by 

r f i — —\ t (- + r (- — - — — 1 r C-a + jl\ 

/^n V 27r»y V2 ~ 2m) \2 2n 2ni J \ 2n ~ 2m) /,,\ 

a2[ ] ~ ~ T ( —) T (- — — + — ) T (- —) T (l — — —) ' 1 ' 

\2ni ) \2 2n^ 2m) \2 2m) \ 2tt 2m j 

This expression coincides with the analogous formula obtained in the GN model for n > 2 
but there is an important difference here: for n < 2, A is negative and therefore the 
amplitudes present no poles in the physical strip. The absence of poles in the physical 
strip matches of course with the nature of massless scattering. It is easy to see that for 
n = 1 (which is simply the massless Ising model) the isoscalar ^-matrix in (|j) correctly 
reduces to —1 whereas for n —>■ 0, the relevant limit for the random bond Ising model, we 
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have 

*""(Mlr(J-A)r(A)) ' ' ' 

In conclusion, our analysis shows that for n < 2 the S-matrix of the elementary fermions 
of the GN model may be regarded as the analytic continuation of the S-matrix for n > 2 
and viceversa. The interpretation, however, is quite different in the two intervals. In fact, 
for n > 2 the S-matrix describes the interaction of massive fermions and, in order to close 
the bootstrap, has to be supported by the computation of the other scattering amplitudes 
relative to the bound states |9|, [10). For n < 2, on the contrary, there are no bound states, 
the fermionic fields split into the chiral components and the S'-matrix describes in this 
case only the interaction of the left-right movers^. 

4. The massless S'-matrix proposed in this paper leaves open some important questions. 
The first question concerns the characterization of the ultraviolet limit of the GN model 
for n < 2. In fact, the perturbative approach to an infrared free theory is usually plagued 
by the presence of Landau pole which may prevent to extend the validity of the theory 
beyond it. If this is the case, the ultraviolet region cannot be reached perturbatively from 
the infrared fixed point. The massless GN S-matrix (|14"D has no singularities for real 
values of 9. This seems to indicate that for the massless GN model, the ultraviolet limit 
is, instead, continuosly connected to the infrared fixed point. In order to clarify this point, 
a first step would be to analyse the behaviour of the massless GN model on an infinite 
cylinder of width R and extract, by means of the Thermodynamical Bethe Ansatz, the 
flow of the effective central charge c(MR) between the ultraviolet (MR — > 0) and the 
infrared (MR — > oo) fixed points. The second question concerns the local properties of 
the GN model and their relationship with the original statistical system described by the 
limit n — >■ 0. As for the limit n — > of the bosonic 0(n) model studied in [14]] in 
connection to the self-avoiding walks, it would be interesting to apply the massless Form 
Factor approach, recently considered in JT5J, and study the behaviour of the correlation 
functions of the random bond Ising model associated to the S-matrix of eq. ([15|) . 
1 Obviously, also the interpretation of the rapidity variable changes in the two intervals. 
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